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ABSTRACT 

We present 10 /im-35 /im Spitzer spectra of the interstellar medium in the Central Molecular Zone (CMZ), 
the central 210 pc x 60 pc of the Galactic center (GC). We present maps of the CMZ in ionic and H2 emission, 
covering a more extensive area than earlier spectroscopic surveys in this region. The radial velocities and 
intensities of ionic lines and H2 suggest that most of the H2 0-0 S(0) emission comes from gas along the 
line-of-sight, as found by previous work. We compare diagnostic line ratios measured in the Spitzer Infrared 
Nearby Galaxies Survey (SINGS) to our data. Previous work shows that forbidden line ratios can distinguish 
star-forming galaxies from LINERs and AGNs. Our GC line ratios agree with star-forming galaxies and not 
with LINERs or AGNs. 

Subject headings: infrared: ISM — ISM: molecules — stars: formation — galaxies: active — galaxies: ISM 
— galaxies: nuclei — galaxies: starburst 



1. INTRODUCTION 

The Galactic center (GC) is the closest gal actic nucleus. At 
a distance of 7.9 ± 0.8 kpc (Re id et alJl2009h . 1 pc in the GC 
corresponds to only 26". The GC provides an opportunity 
to unveil interactions between various physical processes in a 
nuclear environment of a galaxy with excellent spatial resolu- 
tion unapproachable for other galaxies. 

The central 170' x 40' (~ 400 pc x90 pc) region of the 
Galaxy is often called t he Central Molecular Zone (CMZ; 
iMor ris & Serabvrj fl996T) . The CMZ is a massive molec- 
ular cloud complex in the Galaxy, which contains about 
10% of the Galaxy's molecular gas and produces 5%- 
10% of its infrared (IR) and Lyman con t inuum luminosity 
(ISmith et all 1 19781: iNi shimura etafl Il980l: iBally et all Il987[ 
1988; Morris & Serabyn 1996). Earlier radio continuum sur- 
veys revealed that the CMZ is an active star forming region 
and contains the most active star forming cloud in the entire 
Galaxy, the Sagittarius B2 (SgrB2) complex. Herschel obser- 
vations revealed that the dust emission from the C MZ mainly 
arises in a ring-like structure (Molina ri et al.l20lTT) . In this pa- 
per, we use the terms CMZ and GC interchangeably, referring 
to the same ~ 200 pc region in the center of the Galaxy. 

Star formation in the CMZ is inevitably affected by the ex- 
treme physical conditions of the natal clouds, which have an 
order of magnitude higher gas density than in the disk, with 
high gas temperature, pressure, turbulence, strong tidal shear, 
and m illi-Gauss magnetic field strengths (Morris & Serabyn 
1996). Because of these unusual conditions not found in 
nearby normal star-forming regions in the disk, the nature of 
star formation in the CMZ is a subject of active research. 

A key to unlocking the secrets of star formation activities 
in the CMZ is a detailed spectroscopic study of its interstel- 
lar medium (ISM) over a wide range of wavelengths. Mid-IR 
emission lines are particularly useful in analyzing the phys- 
ical properties of the GC gas because it lies behind heavy 
dust obscuration (Ay ~ 30 mag). Forbidden emission lines 
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in the mid-IR are generally insensitive to electron gas temper- 
atures, and can be used to determine physical properties such 
as the electron density and ionization parameters, and to iden- 
tify sources of ionization. Furthermore, molecular hydrogen 
emission from pure rotational transitions can be observed in 
the mid-IR, which may hold vital clues to gas heating mecha- 
nisms in the CMZ. 

A mid-IR spectroscopic survey of the GC was 
previously conducted at 2 /jm-196 pm based 

on ISC; obse rvations ([Rodriguez-Fernandez et al. 

1 200 lal 12004 IRodrfguez-Fernandez & Martin-Pintado 
2005). They observed 15 different lines of sight 
to mo lecular clouds in the CMZ (see Fig ure 1 

in IRodrfguez-Fernandez & Martm-Pintadol 120051) . 

Rodri guez-Fernandez et al.l (1200 lal 120041) used observa- 
tions of molecular hydrogen emission to discuss heating 
mechanisms of warm molecular gas in the CMZ, and con- 
cluded that low-density photon dominated regions (PDRs) 
and low-velocity shocks (v < 10 km s -1 ) are required to ex- 
plain the temperatures derived from the warm molecular gas. 
Rodri guez-Fernandez & Martin-Pin tado (2005) analyzed fine 
structure lines in the GC far from thermal continuum sources 
and massive clusters. They concluded that the ionizing radia- 
tion field is rather constant throughout the CMZ, suggesting 
ionization by relatively hot and distant stars, and found that 
excitation ratios, temperatures, and ionization parameters of 
ionized gas in the CMZ are similar to those found in some 
low-excitation starburst galaxies. 

A more recent survey in the mid-IR was carried out 
by Simpso n" et alJ (120071) using the Infrared Spectrograph 
(IRS; Hou ck et alj|2004l) onboard the Spitzer Space Telescope 
dWerner et al.ll2004 1), whi ch has a higher sensitivity than ISO 
at 10 /im-38 /im. The Si mpson et al.l survey consists of 
high-resolution spectra of 38 positions along a narrow 24' 
long strip at a Galactic longitude / « +0.1°. iSimpson et alJ 
measured several forbidden emission lines and molecular hy- 
drogen lines to constrain the physical conditions of clouds 
near the Quintuplet Cluster, the Arches C luster, the Radio Arc 
Bubble (Rodriguez-Fernandez et al. 2001), and Arched Fila- 
ments. They concluded from their observations that the main 
source of excitation in the GC is photo-ionization from the 
massive star clusters and that multi-component PDR models 
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can explain th e observed line emission from m olecular hydro- 
gen (see also IContinil 120091) . ISimpson et al.l also concluded 
that shocks in the Radio Arc Bubble are responsible for strong 
[O IV] emission. 

In our recent spectroscopic surv ey of massive youn g stel- 
lar objects (YSOs) in the CMZ (lAnet al.l 12009. 201 3), we 
used Spitzer/IRS to collect an extensive set of mid-IR spec- 
tra for 107 YSO candidates. The goal of this survey was to 
discover and characterize the spectroscopic properties of mas- 
sive YSOs in the GC. To achieve this original goal, we spent 
half of our observing time on background spectra near each 
YSO candidate because of the strong and spatially variable 
background emission in the GC. Our GC background spec- 
tra, which are the by-product of our YSO observing program, 
now constitute the largest and most comprehensive mid-IR 
spectroscopic data set available to study the properties of the 
ionized and molecular gas in the star-forming nucleus of the 
Galaxy. 

In this paper, we present panoramic mapping results from 
ionic and molecular hydrogen emission lines throughout the 
entire CMZ. Data acquisition and spectral analysis are pre- 
sented in §|2] Line intensities and radial velocities in the CMZ 
are presented in § [3] Mid-IR line ratio diagnostics are used 
to compare the physical properties of the Galactic nucleus to 
those of other nearby galaxies in § [4] Our results are summa- 
rized in § |5] 

2. METHOD 

2. 1 . Observations and Extraction of Spectra 

The IRS spectra presented in this paper were obtained in 
May and October 2008 as part of Spitzer Cycle 4, during 
our observing program to identify massive YSOs in the GC 
(Program ID: 40230, PI: S. Ramirez). We targeted 107 point 
sources with extremely red colors on nea r- and mid-IR co lor- 
color diagrams as YSO candidates (see lAn et al.l 1201 Tl for 
details). The positions of our spectra are shown in the top 
panel of Figure Q] on top of the Spitzer /IRAC 8 ^m image 
dStolovv et ai1l2006t iRamfrez et ail 120081) . Our survey en- 
compasses a large area in the flattened CMZ cloud complex 
(\b\ < 0.24°), covering regions near strong radio continuum 
sources (Sgr A, Sgr Bl, Sgr B2, and Sgr C), massive star 
clusters (the Quintuplet, Arches, and Central clusters), and 
the Radio Bubble (see the bottom panel of Figure[T). 

In our Spitzer program, massive YSO candidates were ob- 
served using both high- and low-resolution IRS modules. In 
this work, however, we utilize only the high-resolution ob- 
servations taken with the short-high (SH; 9.9 /im-19.6 pm, 
A/AA-600;4.7" x 11.3" or 0.18 pc x0.43 pc slit entrance) 
and the long-high (LH; 18.7 /xm-37.2 pm, A/AA ~ 600; 
11.1" x 22.3" or 0.43 pc x0.85 pc slit entrance) modules. 
These two modules share the same slit centers on the sky. 

Four background spectra were taken at ~ 1' (2.3 pc) away 
from each of the YSO candidates; the positions of these points 
can be glimpsed by looking at the open circles in the top 
panel of Figure Q] each of which has a 1' radius. The ex- 
act locations of these spectra were carefully chosen by visual 
inspect ion to avoid bright point sources on the IRAC 8 pm 
image (Sto lovv et alJ2006l ; IRamfrez et al.l2008l) . and were in- 
tended to provide a level of background mid-IR emission sim- 
ilar to that of the target position. As shown in Figure Q] our 
background spectra are somewhat uniformly distributed over 
a - 90' x 25' (~ 210 pc x60 pc) region in the CMZ. 

We reduced IRS spectra from the basic calibrated data 



(BCD) products version SI 8 . 7 . 0. We corrected for rogue 
pixel values using the software package IRSCLEA1S0 pro- 
vided by Spitzer Science Center (SSC), and applied the 
DARKSETTLE software package to the LH frames to correct 
for non-uniform dark current. We extracted SH and LH spec- 
tra using the SPICE tool in an "extended" extraction mode 
(with slit loss corrections for a source infinite in extent), and 
further corrected for fringe patterns using the IRSFRINGE 
package. More information on our IRS obse rvations and ba- 
sic data reduction is found in I An et aTl (1201 ll) . 

The contribution of zodiacal light is relatively high in the 
GC, typically amounting to ~ 10%-30% of the continuum 
emission from the ISM in the IRS spectral range. We used the 
zodiacal light estimator in the SSC-provided software SPOT, 
which is based on COBE/DIRBE measurements, to determine 
the contribution from the zodiacal light in each season. The 
zodiacal spectrum has a peak emission of ~ 40 MJy sr" 1 at 
20 pm, with a ~ 10% seasonal variation. We constructed 
a smoothed zodiacal spectrum using quadratic interpolation, 
and then subtracted it from the extracted spectra. Since our 
IRS targets are found within a degree of the GC, spatial vari- 
ations in the zodiacal emission are negligible. The correction 
for zodiacal emission has no direct impact on emission line 
flux measurements, but has an influence on the extinction cor- 
rection since we measure foreground extinction from t he co n- 
tinuum emission near the 10 pm silicate feature (see § !2.31 l. 

We used the above procedure to obtain ISM spectra from 
SH and LH for 428 individual lines of sight in the GC (i.e., 
four background spectra for each of 107 point-source tar- 
gets). Figure [2] shows one of the observed spectra in the 
GC, after the zodiacal light correction. Different IRS spec- 
tral orders are shown in alternating colors. Forbidden emis- 
sion lines and molecular hydrogen lines are marked, and the 
wavelength ranges of strong poly cyclic aromatic hydrocarbon 
(PAH) emission features are indicated. 

We identified spectral orders that best display individual 
lines, and used Gaussian profile fits to compute line fluxes, 
as described in the following section. Figure [2] show s that or- 
der til ts are present in some of the IRS spectra. In lAn et al.l 
d201 ll) . we used low-resolution spectra to correct for order 
tilts before we combined individual high-resolution spectra 
from different orders. In the following analysis, however, we 
combined spectra without order-tilt corrections, because our 
low -resolution background slits are not always co-spatial with 
the high-resolution slit positions. Figure [2] shows that the line 
fluxes should not be strongly affected by the order tilts. 

2.2. Measurement of Emission Line Strengths and Radial 
Velocities 

We measured emission line fluxes from several ionic fine- 
structure lines, as well as molecular hydrogen lines from pure 
rotational transitions, H 2 0-0 S(0), H 2 0-0 S(l), and H 2 0-0 
S(2). A list of emission lines included in this paper is shown 
in Table Q] in order of increasing wavelength. In addition to 
these lines, we attempted to measure a flux from other ionic 
features such as [Ar V] 13.1 pm, [P III] 17.89 pm, [Fell] 
17.94 pm, [Fe II] 24.52 ^m, and [Fe II] 35.35 pm, but none 
of these features were strong enough to be detected. The third 
column in Table Q] shows IRS modules and spectral orders, 
from which individual line fluxes were measured (see below). 
The fourth column shows normalized extinction coefficients 

4 The SSC software packages can be found at 

|http: //irsa. ipac ■ caltech . edu/data/SPITZER/| 
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FIG. 1 . — Top: Spatial distribution of 107 SpitzerflRS targets on the IRAC 8.0 pm image I Stolovy et al. 2006; Ramirez et al. 2008) of the Galactic center (GC). 
The above image sh ows the entire Centra l Molecular Zone (CMZ), which covers approximately 170' X 40' (~ 400 pc x90 pc) centered on the GC. Our IRS 
point-source targets (An et al. 2009, 201 1) are marked with 1' (2.3 pc) radius circles. In this paper, we analyze the GC interstellar medium (ISM) as measured by 
four background spectra that are ~ 1' away from each target. Bottom: Schematic diagram of locations of the CMZ molecular complexes (Sgr A, Sgr Bl, Sgr B2, 
and Sgr C; large open circles), the Radio Bubble (large open circle), and star clusters (Quintuplet, Arches, and Central clusters; open star symbols). The grey 
circles represent the location of the Spitzer/IRS targets. 
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from the GC extinction curve in lChiar&Tielensl(l2006b : see 
also the bottom panel of Figure|2] The ionization potential for 
each ion is listed in the last column in TableQ] 

Figure|3]shows examples of the line profile fits for the emis- 
sion lines listed in Table Q] Grey histograms are observed 
IRS spectra from various lines of sight. The red line repre- 
sents the best-fitting Gaussian profile, obtaine d using the non - 
linear least squares fitting routine MPFIT (lMarkwardtl l2009 ) . 
The underlying blue line shows a local continuum constructed 
from a l "-3' r/ order polynomial fit to the continuum points on 
each side of the emission line. A standard deviation of data 
points from the continuum line was used as an effective lcr 
error per data point in the profile fitting. High order poly- 
nomials were used to determine a local continuum near PAH 
emission features. Simultaneous fits of two Gaussian profiles 
were made for two sets of blended emission lines, [Ne V] 
14.32 fim and [CI II] 14.37 ^im, and [O IV] 25.89 ^m and 
[Fell] 25.99 fim. 

The total line flux was computed by integrating the under- 
lying flux of the best fitting Gaussian profile in wavenumber 
space. Errors in these line fluxes were estimated by adding in 



quadrature the 1 a flux uncertainties derived from uncertain- 
ties in the height and width of the Gaussian fit, or by propa- 
gating lcr errors in the local continuum, whichever is larger. 
The latter was computed as 

<7flux ~ ^rms A^ max ~ ^Vms ^Tl J ^^ max ' 

where F nns is the rms dispersion in the local continuum. The 
parameter AA max is the maximum value of the full-width at 
half-maximum (FWHM) of the line profile set in the Gaussian 
line fitting, which we assumed to be Av = 600 km s" 1 for 
all lines. In most cases, errors propagated from Gaussian fits 
were larger than those estimated using EquationQ] 

The line center was allowed to vary in our line fitting pro- 
cedure, from which we measured line-of-sight velocities (v r ) 
from individual emission lines. Although high-resolution 
spectra from Spitzer/IRS have a relatively low spectral res- 
olution (R = A A/ A = 600), radial velocities can be obtained 
with a precision of a few t ens of kilometers per second for 
strong emission lines (e.g.. ISimpson et all 120071) . We com- 
puted radial velocities in the Local Standard of Rest (LSR), af- 
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I Observed ISM spectra (l,b) = (+0. 259, -0.012) \ 
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FIG. 2. — Top: Observed SpitzerflRS spectrum in the ISM of the GC, measured at (Lb) = (+0.2588°,— 0.0119°), after subtracting the zodiacal emission. 
Different spectral orders are shown in alternating colors. Ionic forbidden emission lines and molecular hydrogen lines are marked with vertical lines, with dotted 
lines indicating non-detected emission lines. Vertical grey strips centered at 10.2 fim and 13.9 fim mark continuum regions used to estimate th e 9.7 fim silicate 
optica l depth (rg 7) for each line of sight (see text). Bottom: Extinction curve for the GC, normalized to the extinction value in the K passband (Chiar & Helens! 
2006). 



TABLE 1 

Emission Lines Detected in Spitzer/IRS Spectra 



Line 


Wavelength 


IRS Module/Order" 


A x /A K b 


IP C 




(jim) 






(eV) 


[S IV] 


10.51 


SH/19 


1.080 


34.79 


H 2 S(2) 


12.28 


SH/17 


0.478 




HI 7-6 


12.37 


SH/17 


0.462 




[Ne II] 


12.81 


SH/16 


0.435 


21.56 


[Ne V] 


14.32 


SH/14 


0.430 


97.12 


[CI II] 


14.37 


SH/I4 


0.431 


12.97 


[Ne III] 


15.56 


SH/I3 


0.533 


40.96 


H 2 S(l) 


17.04 


SH/12 


0.698 




[S III] 


18.71 


SH/I1 


0.762 


23.34 


[Fe III] 


22.93 


LH/17 


0.622 


16.19 


[Ne V] 


24.32 


LH/16 


0.595 


97.12 


[O IV] 


25.89 


LH/15 


0.562 


54.94 


[Fe II] 


25.99 


LH/15 


0.560 


7.90 


H 2 S(0) 


28.22 


LH/14 


0.513 




[S III] 


33.48 


LH/12 


0.433 


23.34 


[Si II] 


34.82 


LH/11 


0.426 


8.15 



a IRS modules and spectral orders used in the extraction of line fluxes. 

Extinction curve in IChiar & Tielemil<200eil) . 
c Ionization potential for each ion. 

ter correcting for the spacecraft motion: the correction terms 
(keyword VLSRCORR provided by the SSC) for May and Oc- 
tober 2008 observing runs were +38.9 km/s and -18.5 km/s, 



respectively. 

The accuracy of our radial velocity measurements was eval- 
uated as follows. First of all, we compared v r measurements 
from the same ionic species, [S III] 18.71 and [S III] 
33.48 fim, which are found on two different IRS modules. 
We found that the mean difference in v r between these lines 
is Av r = 31 km s , indicating the error from a module-to- 
module and/or order-to-order change in the IRS wavelength 
calibration. In addition, the scatter in the v r difference from 
these two lines is a v = 17 km s" 1 , which is a measure of the 
precision in the v r determination. 

The size of systematic errors in v T can be further examined 
by comparing v, determined among strong ionic emission 
lines. We used [Nell] 12.81 fim as a basis of comparison in v r . 
We found a median difference in v r for all of the IRS spectra at 
I < +0.2° as follows: Av r = +5.3 km s" 1 for [Ne III] 15.56 fim, 
Av r = -10.8 km s" 1 for [S III] 18.71 fim, Av r = +21.7 km s" 1 
for [S III] 33.48 fim, and Av r = +71.1 km s" 1 for [Si II] 
34.82 fim, respectively, in the sense that a positive differ- 
ence indicates a larger v r from a given line than from [Ne II]. 
We used these Av r values to first shift all ionic lines to the 
[Ne II] velocity. The [Ne II] line has on average a smaller v r 
by Av r = 22.60 km s" 1 than the mean vlsr = +51.0 km s" 1 
measured from a hy drogen recombination line in Sgr Bl 
dMehri nger et al] | 19921) . We thus a pplied this second corr ec- 
tion to shift all ionic lines to match Mehring er et all (119921) in 
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FIG. 3. — Examples of line profile fits for the GC ISM spectra. Grey histograms show observed Spitzer/IRS spectra, displayed over the same wavelength 
interval, but with an arbitrary intensity scale on the vertical axis. Ionic emission lines as well as pure rotational lines from molecular hydrogen, such as H2 — 
S(0), H2 0-0 S(l), and H2 — S(2), are shown. The best-fitting Gaussian profile is shown as a red line on top of a local continuum (blue line). 



SgrBl. 

Molecular hydrogen lines from pure rotational transitions 
are also strong enough to be detected in the GC and provide 
radial velocities. We compared v r from individu al IRS spec- 
tra to the peak v> of CO J = 4 — > 3 emission (Mart in et al.l 
120041) at -0.2° < I < 1.0°, and found Av r = +62.5 km s _1 , 
+51.2 km s" 1 , and +29.4 km s" 1 , for H 2 S(2) 12.28 /im, H 2 
S(l) 17.04 /xm, and H 2 S(0) 28.22 /mi, respectively. The 
sense of the difference is that radial velocities from these lines 
are on average lower than those from the CO line. We then 
shift all H 2 lines to the CO velocity. 

The above systematic offsets in the measured v r indicate 
that the wavelength calibration error in the Spitzer/IRS spec- 
tra is on the order of at least a few tens of kilometers per sec- 
ond. To summarize, we put all v r measurements from the IRS 
observations in the Local Standard of Rest (LSR) by applying 
zero-point offsets to vy. 



2.3. Foreground Extinction Estimates 

We corrected line flux measurements for foreground ex- 
tinction between the GC and the Sun. This was done on a 
spectrum-to-spectrum basis, because of the patchy dust ex- 
tinction in the GC. We utilized two different approaches to 
correcting for foreground extinction. 

The first approach is based on the method developed in 
Simps orTet al.l ([2007), who determined the flux ratio at 10 /mi 
and 14 /im {Fu/Fiq) and used it to infer the optical depth 
of the 9.7 /im silicate absorptio n feature (rqj) near the Ra- 
dio Arc Bubble. Following the ISimpson et al.l prescription, 
we estimated mean fluxes at 10.00 /im < A < 10.48 /im and 
13.50 /im < A < 14.30 /im from SH spectra after a 3a re- 
jection. As shown by vertical grey strips in Figure |2] these 
two continuum wavelength ranges are almost free of PAHs 
features and ionic emission lines. Then we estimated T9.7 for 
each line of sight, by using a simple linear relationship be- 
tween the 9.7 /im silicate optical depth and the continuum 
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flux ratio, 7-9.7 = (ln(Fi 4 /F 10 ) - 0.809)/0.560 which was in- 
ferred from tabulated values in Simpson et al. (2007). The top 
panel in Figure|4]shows the spatial distribution of 79 .7 derived 
in this way for all lines of sight in our program. 

Al ternatively, we utilized the A v extinction map toward the 
GC (ISchultheis et al. 2009) based on the observed stellar lo- 
cus of red giant stars on 2MASS and SpitzerflRAC color-color 
diagrams. T he bottom panel i n Figu re|4]shows the distribution 
of 79 .7 from Schulth eis et al.l (120091). after converting thei r Ay 
values into 79.7 using Ay / Tgj = 9 (Roche & Aitk enlll985l). 

Figu re [5] compares 79 7 derived using the Simpson et al. 
(2007) method to the visual extinction (Ay) from 
Schul theis et al.l ( l2009t) . Figure [5] shows a large scatter 
between these two independent extinction estimates. How- 
ever, their mean trend agrees well with the GC relation 
(lRoche& Aitkenl[T98ll Ay/rg? = 9). For comparison, we 
also s how the relationship for the local ISM (Roc he & Aitkenl 
11984 Ay/79.7 = 18.5) in Figure [5] which does not match the 
observed trend in the GC at all. 

In the fol l owing analysis, we derived 797 using 
Simps orTet al.l (120071) method and present results after 
correcting line fluxes for extinction, unless otherwise stated. 
However, we r epeated our ana l ysis w ith foreground extinc- 
tion from the ISchultheis et al.1 (120091) map, and compared 
results with each other. As shown below, our main results are 
solid, and are only weakly dependent on the specific choice 
of foreground extinction correction. 

In addition to the above extinction corrections, we further 
corrected observed line fluxes from [Ne III] 15.56 /im for 
the absorption induced by CO2 ice grains on a spe ctrum-to- 
spectru m basis. We used the same procedure as in lAn et al.l 
d201 ll) to decompose a wide CO2 ice absorption band over 
~ 15 pm-16 /im using five laboratory spectral components 
of different CO2 ice mixtures, and obtained an optical depth 
of the CO2 ice from the best-fitting model profile. The mean 
optical depth from all of the GC ISM spectra is (At) = 0.065 
at the position of the [Ne III] line. The relatively small opti- 
cal depth at 15.56 /im is because the [Ne III] line is located 
in the long-wavelength wing of the CO2 ice absorption, in ad- 
dition to the fact that the CO2 ice band in the ISM i s weaker 
and n arrower than those seen in massive YSOs (see lAn et al.l 
l20Tlh . 

For a given 79.7, we corrected a line flux for interstellar ex- 
tinction at each wavelength using 

where A\/Ak is the GC extinction curve value in 

IChiar & Tielensl (120061) . normalized to the extinction in the 
K bandpass (see Table [TJ. We adopted Ak/Av rs 0.11 
dFiger et al.|l 19991) and the GC rela tion of A v /r 9 . 7 = 9 as deter- 
mined by Roche & Aitken ( 1985). The extinction-corrected 
line flux was then obtained using 

log/o = log/ obs + 0.4A A . (3) 

2.4. Coadded GC ISM Spectra 

In addition to individual ISM spectra, we also created and 
analyzed a coadded GC spectrum as shown in Figure [6] This 
allows us to compare the ISM spectrum measured across 
~ 200 pc to the individual spectra with the ^ 3.5 pc resolu- 
tion. We constructed the coadded IRS spectrum by summing 
fluxes from 428 individual high-resolution spectra in the GC. 



The top panel in Figure [6] shows the spectrum coadded with- 
out any foreground extinction correction. The bottom panel 
displays the result of c orrecting each individ ual spectru m fo r 
extinction based on the lSimpson et al.l (12007) method (§ !2.31 l. 
then coadding the corrected spectra. We also illustrate the 
spectrum created by correcting each i ndivid ual spectrum for 
extinction with the Schulthei s et al.l d2009) extinction map, 
then coadding. The coadded spectra from each extinction cor- 
rection method are very similar. We present below line fluxes 
from the coadded spectrum together with those obtained from 
individual spectra. 

3. RESULTS 
3.1. Panoramic Emission Line Mapping in the GC 

Panoramic emission line maps are displayed in Figure |7]in 
a logarithmic flux scale for each line. Only those lines de- 
tected at more than a 3cr level were included. Table|2]provides 
the individual line fluxes us ed in these maps, after correcting 
for dust extinction using the Simpson et al. (2007) technique. 
Only a portion of Table [2] is shown here to demonstrate its 
form and content, and a machine-readable version of the full 
table is available in the electronic edition of the Journal. 

Table [3] lists the line fluxes measured by coadding all GC 
spectra, with three different extinction correction methods. 
The line fluxes and ratios listed in the second column are 
those measured from a spectrum coadded from individual 
spectra that ha v e firs t been corrected for extinction using the 
Simps on et al.l (120071) F14/F10 technique. Similarly, values 
in the third column are those measured in a spectrum coad- 
ded from individual s pectra that were first corrected using the 
Schul theis et al.l (120091) extinction map. The last column lists 
line fluxes and ratios measured by coadding all GC spectra, 
then applying an extinction c orrection for T9 . 7 = 3. 439 (the 
median value obtained from the Simpso n et al.l (120071) Fn /Fm 
technique; see Figure|5]l. 

Mapping results in Figure |7]are shown in the same order of 
increasing wavelength as in Table [TJ Each pixel in the map 
covers a 1.5' x 1.5' region of the sky (~ 3.5 pc x3.5 pc), 
which is an order of magnitude larger than the area covered 
by a slit entrance of either SH or LH. We divided each ex- 
tracted line flux by the areal coverage of the corresponding 
slit entrance (53 arcsec 2 for SH and 248 arcsec 2 for LH), and 
computed unweighted mean intensities (W m" 2 sr" 1 ) in each 
pixel of the map in Figure [7] Each pixel includes two back- 
ground IRS pointings on average. 

Average line intensities in Figure [7] we re corrected for fore- 
ground du st ext inction derived using the Sim pson et al.l (120071) 
method (§ I2.3t . which relies on the continuum flux ratio be- 
tween 10 pro. and 14 pm (see the top panel in Figure 0J. 
Mid-IR forbidden emission lines, such as [Ne II] 12.81 pm, 
[Ne III] 15.56 pm, [S III] 18.71 /mi, [S III] 33.48 /mi, and 
[Si II] 34.82 /im, are strong in the CMZ. Their line strengths 
vary across the region, with strongest emission observed near 
the Arches cluster and Sgr Bl. Weaker emission is observed 
in the CMZ for H I 7-6 12.37 pm, [CI II] 14.37 /mi, [Fe III] 
22.93 /im, and [Fe II] 25.99 pm. Emission from pure rota- 
tional H 2 lines, 0-0 S(0), 0-0 S(l), and 0-0 S(2) at 28.22 pm, 
17.04 /im, and 12.28 /im, respectively, is also strong, but these 
H2 lines are relatively constant ove r the CMZ. We will dis- 
cuss the uniform H2 emission in § 13.21 further utilizing ra- 
dial velocities measured from the IRS spectra. Fine structure 
lines from highly ionized species such as [O IV] 25.89 /im 
are found throughout the CMZ, but [S IV] 10.51 p,m, [Ne V] 
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TABLE 2 

Line Flux Measurements from Individual Spectra 



Galactic 


Galactic 


[Ne 11] 


H 2 S(l) 




[S III] 


H 2 S(0) 


[S III] 




[Si II] 


longitude 


latitude 


12.81 fim 


17.04 fim 


1 


8.71 ^m 


28.22 fim 


33.48 fim 


34.82 ^m 


(deg) 


(deg) 


















-0.7954 


-0.0263 


7.32 ±0.11 


8.67 ±0.08 


2 


87 ±0.15 


1.99±0.15 


4.25±0.14 


16 


88 ±0.20 


-0.8254 


-0.0419 


10.94±0.13 


8.10±0.09 


10 


21 ±0.17 


1.87±0.13 


16.04 ±0.35 


25 


87 ±0.45 


-0.8037 


-0.0450 


4.99 ±0.11 


6.83 ±0.08 


2 


01 ±0.15 


1.68 ±0.05 


3.42 ±0.11 


11 


57±0.16 


-0.8198 


-0.0212 


8.10±0.11 


6.65 ±0.07 


3 


99 ±0.22 


1.46±0.10 


6.64 ±0.11 


15 


18 ±0.32 


-0.7488 


-0.0254 


4.43 ±0.13 


5.13±0.18 


1 


91 ±0.23 


1.60±0.10 


3.68±0.18 


7 


82±0.16 


-0.7814 


-0.0477 


5.51 ±0.20 


9.42 ± 0.42 


2 


93 ±0.34 


2.08 ±0.07 


3.54±0.31 


11 


85 ±0.19 


-0.7594 


-0.0510 


6.50±0.18 


9.64 ±0.56 


3 


41 ±0.54 


2.75 ±0.15 


4.65 ±0.40 


13 


44 ±0.27 


-0.7720 


-0.0227 


7.00 ±0.1 8 


10.87 ±0.33 





00 ±0.00 


2.38 ±0.09 


5.50 ±0.34 


14 


19 ±0.21 


-0.7529 


-0.0297 


25.84 ±0.70 


54.67 ±0.95 


44 


04 ±2.32 


13.19±0.60 


28.87 ±1.07 


58 


81 ±1.89 


-0.7836 


-0.0515 


6.07 ±0.10 


9.02 ±0.09 


2 


38±0.15 


1.98 ±0.07 


3.50 ±0.59 


12 


99 ±0.71 



NOTE. — Table 2 is published in its entirety in the electronic edition of the Asrrophysical Journal Supplement Series. A portion is 
shown here for guidance regarding its form and content. The electronic edition of Table 2 contains the observed fluxes and errors for the 
emission lines listed in Table 1. 



14.32 /im and [Ne V] 24.32 fim were detected only in a few 
lines of sight to the GC. Our mapping resu lts and interpre- 
tation remain qualitatively unchanged if the Schulthei s et alj 
(2009) map is used for the foreground extinction correction. 



3.2. Molecular Hydrogen Line Emission and Radial 
Velocity Mapping 



Figure [7] shows line intensity maps for pure rotational tran- 
sitions from the lowest three levels of molecular hydrogen, 
H 2 S(2) 12.28 fim, S(l) 17.04 /xm, and S(0) 28.22 /mi. We 
detected these lines in almost all lines of sight in the GC at 
more than a 3er level. Their intensity distributions are rather 
uniform in the CMZ, compared to the spatial structures seen 
for ionic forbidden emission lines. We also find weak correla- 
tions in line intensities between the pure rotational hydrogen 



8 



An et al. 



TABLE 3 

Line Flux Measurements from Coadded Spectra 







Extinction Correction 




Line Flux 


Simpson et al. (2007) 11 


Schultheis et al. (2009)" 


Posterior Correction 


(Line Ratio) 


(lir 7 Wm -2 sr 1 ) 


(10- 7 Wm- 2 sr 1 ) 


(icr 7 Wm" 2 sr 1 ) 


1VJ 1U.51 (JXW 


2. 30 ± 3.22 


O.UO it 2.yu 


U.54 ± 4.38 


H2 o{l) 12. zo fim 


/.06 ± 0.34 


8.45 ± 0.34 


£. iq _i_ 1 *;n 
6.28 ± 1.60 


ITt"7 £. 1T n . 

H I 7-6 12.37 fim 


1 HI 1 f\ 1 "7 

1.67 ± 0.17 


1.95 ±0.17 


1.36± 0. /5 


nvT„ ttT 11 01 . 

[JNe IIJ 12.81 fim 


CO (\C\ 1 A TA 

58. oy ± O.JO 


/if Art I n ti 

63.00 ±0.32 


lie n<i i 1 tj 

45. y4 ± 1.23 


[NeV] 14.32 fim 


0.00 ±0.08 


0.00 ±0.08 


0.06 ±0.32 


rpi ir"| 1 4. 37 iim 
[V^l 11J Jt. J / fllil 


n an _i_ n on 
U.JZ _n u.zu 


n it. _|_ n OQ 

U.JJ __L \J.\) J 


22 ± 37 


[Ne III] 15.56 Aim 


4.74 ±0.26 


4.94 ±0.26 


3.80±1.48 


ri2 I J I / .U'f fim 


2U. J / it U.2j 


ia on 4- n in 


1 ^ Al 4- 1 IS 
1 It Z.Zo 


[3 IIIJ to. / 1 /iDl 




yy. ji ± u.o / 


jz.yi it z.j4 


rc a tit"! TO AT ,, m 

[re IIIJ zz.yj £tm 


/i 1 1 _i_ n 1 q 
4. 12 it 0. Ly 


4.UJ ± U. 15 


z.oy it i.jy 


[INe VJ Z4.jz /iDl 


U.23 ± o. iy 


mixing 


U. 10 it 1 .oU 


[0 IV] 25.89 fim 


0.70±0.12 


0.86±0.13 


0.60 ±0.45 


[Fe II] 25.99 fim 


2.97 ±0.12 


3.48±0.14 


2.26 ±0.83 


H 2 S(0) 28.22 fim 


3.25 ±0.59 


3.85 ±0.61 


2.73 ± 1.27 


[S III] 33.48 /xm 


75.30 ±0.98 


77.48 ±1.03 


58.07 ±0.77 


[Si II] 34.82 fim 


61.76 ±1.34 


70.83 ±1.41 


51.29±1.04 


[Ne III] / [Ne II] 


0.082 ±0.005 


0.078 ± 0.004 


0.083 ±0.032 


[Si II] / [S III] 


0.820 ±0.021 


0.914±0.022 


0.883 ±0.021 


[Fe II] / [Ne II] 


0.051 ±0.002 


0.055 ±0.002 


0.049 ±0.018 


[0 IV] / [Ne II] 


0.012 ±0.002 


0.014 ±0.002 


0.013 ±0.010 


[Fe II] / [O IV] 


4.228 ±0.766 


4.071 ±0.624 


3.766± 3.120 



a Fluxes corrected for dust extinction using the Simpson et al. i 2007) technique (see text). 

Fluxes corrected for dust extinction using the Schultheis et al. ( 2009) map (see text). 
c Extinction corrections applied to the coadded spectrum assuming {T9.7) = 3.439 (see text). 



a combination of PDRs and diffuse ionized gas can be used 
to explain the obser ved dust, H?, ne utral gas, and ionized gas 
emission in the GC. lPak et alJ(l!996h also found that the most 
likely cause of the large-scale ro-vibrational H2 v= 1 — !>0S(1) 
emission at 2.12 /xm in the GC is UV excitation by hot mas- 
sive st ars. However, Rodriguez-Fernandez & Martin-Pintado 
(120051) showed that while PDR models can reproduce the ob- 
served pure rotational H2 lines from excited levels in the GC, 
they are not enough to explain all the emission from the lowest 
levels, S(0) and S(l). They argued that low velocity shocks or 
turbulent motions are needed as an additional heating mech- 
anism to reproduce excess emission from low e xcitation H2 
lines. On the other hand, Simps orTet all d2007) concluded, 
based on their H2 rotational line measurements, that multi- 
component models of warm molecular gas in PDRs along the 
line of sight to the GC can fully explain the observed H2 line 
ratios, without requiring shocks. 

To further investigate the source of low-level rotational H2 
lines, we compared radial velocity distributions of various 
mid-IR emission lines. Figure [8] shows panoramic radial ve- 
locity (v r ) maps of the CMZ for several strong ionic and 
molecular hydrogen emission lines. Radial velocities from 
strong ionic lines exhibit a systematic rotation of the CMZ, 
consistent with the rotation of the Galactic disk. The sense of 
the rotation is that the eastern part of the GC, including the 
Sgr B complex, is systematically receding from the Sun (pos- 
itive v r ), and the Sgr C complex is systematically approaching 
the Sun (negative v r ), with respect to the dynamical center of 
the Galaxy, indicated by a vertical dotted line at Sgr A*. 

Radial velocity distributions as a function of Galactic lon- 
gitude are displa yed in Figure [9] on top o f the l-v r diagram 
of CO 7 = 4—^3 dMartin et al.ll2004l) : see § EZlfor details on 
the v r determination. Only those v r measurements near the 
Galatic plane, \b\ < 0.1°, are included in the l-v r diagrams, 



05 
O 
O 

C\} 

- 60 - 

CO 



m 40 

CD 



3 

o 

CO 



20 



< 



1 1 1 1 1 1 1 

Local/ ISM 


1 | , 1 1 | 

/ 


: / " 


Galactic 
center 




fe" xx ~ 


*y 'vs^^^xxx* x 


* xjj 

J5 

XX X 


/ * /'x* * ** " 

" ,,1,1,1, 


, 1 , , , f 



2 

T 



9.7 



4 6 
(Fu/F 10 ) 



8 



n pai 

derived using the Simpson et al. ( 2007) method (F14/F10) and the visual ex- 
tinction (Ay) from the Schultheis et al. (2009) extinction map for all of the 
IRS background spectra. Dashed lines are previously me asured relationships 
in the line of sight to the G C (A v /t 91 = 9;|R oche & Aitken 1985) and for the 
local ISM (A v /t 9 j = 18.5; Roche & Aitken 198j), respectively. 



emission and ionic forbidden lines (0.3 < p < 0.6); correla- 
tions are stronger (p > 0.9) among strong forbidden emission 
lines. 

The pure rotational H2 lines from warm molecular gas are 
importan t tools for studying heating mechanisms in the GC 
(e.g., see Rous sel et al.ll2007[, for extragalacti c H2 line mea- 
surement). Rodriguez-Fernandez et al. (2004) suggested that 
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FIG. 6. — Coadded GC ISM spectra constructed from 428 high-resolution IRS spectra. Locations of ionic forbidden emission lines and molecular hydrogen 
lines are marked, where dotted lines indicate positions of weakly detected emission lines. Top: Coadded spectrum created by adding fluxes from individual 
IRS s pectra without foregro und extinction corrections. Bottom: Same as in the top panel, but coadded after correctin g individual spectra for extinction based 
on the|Simpson et al. 1 2007) method. Underlying grey line shows the same but with extinction corrections based on the Schultheis et al. (2009) extinction map. 
Different spectral orders are shown in alternating colors. 



and averaged antenna temperatures (7^*) from the CO survey 
are shown in the same latitude range. 

Three large cross signs in the upper panel of Figure |9]in- 
dicate radial velocit i es from radio recombination-line studie s 
ddePree etal.lll996t lLang et al.ll200lt iLiszt & Spikerl jl995). 
independent of the v r measurement in Mehringer e t al.l(11992l) . 
These show an excellent agreement with our v r values from 
ionic lines (middle panel). The IRS ionic lines generally fol- 
low a lower branch (smaller v r ) in the CO l-v r diagram, and 
the highly negative radial velocities observed in the northern 
rim of Sgr A are associ ated with the Arche d Filaments and 
the Radio Bubble (e.g., lSimpson et al.ll2007h . The IRS spec- 
tral resolution is not high enough to resolve individual radial 
velocity structures along the line of sight as in the CO sur- 
vey, and our v r measurements refer to those that belong to the 
highest peak of the emission. 

Radial velocities from the molecular hydrogen lines S(0), 
S(l), and S(2) are displayed in the bottom panel of Figure [9] 
As shown in this panel, H2 lines show a systematically flat- 



ter v r curve than those from forbidden lines (middle panel). 
The flatness of the v r curve depends on H2 excitation, with 
S(0) showing the most difference from [Ne II] and S(2) show- 
ing the least. The H2 velocities and [Ne II] velocities differ 
most in the region betw een the Arches cluster and Sgr A, 
and toward Sgr C. Simpson et al. (2007) also found a sys- 
tematically different v r between ionized species and molec- 
ular hydrogen. Our mapping results show that v r from the 
lowest H2 energy level, S(0), even exhibits flat rotation near 
Sgr A*. A flat v r distribution in the longitude vs. radial ve- 
locity diagram toward the GC is typically attributed to gas 
clouds in the disk along the line of sight to the GC (e.g., 
iBinnev et alj|1991t iRodrfguez-Fernandez et all 120061) . This 
suggests that the H2 line-emitting clouds are a superposition 
of several dense clouds along the line of sight to the GC. The 
v r offsets for H2 emission lines are still un certain, because we 
opted to use the CO J = 4 -+ 3 data cube dMartin et al.ll2004l) 
at a particular longitude range, -0.2° < I < +1.0° (see § 12.21 ). 
Both the uniform H2 intensity distributions and distinct H2 
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Figure [T}. 



radial velocity structures suggest that a significant fraction of 
the S(0) emission, and a smaller fraction of S(l) and S(2) 
emission, arise in clouds that are probably not associated with 
the mid-IR forbidden line emitting cloud complexes in the 
CMZ. 

3.3. Emission Line Ratio Mapping in the GC 

Mapping results for forbidden line ratios are presented in 
Figure [T0| Line flux es with extinction corrections from the 
Simpson et al. (2007) Fu/Fw technique were used in the 
computation of these line ratios, and only those detected at 
more than a 3cr level were included in Figure [10] A moving 
boxcar average for these line ratios is shown on the right pan- 
els as a function of the Galactic longitude, where the error bars 
indicate the standard deviation of data points in each longitude 
bin. We also present moving bo xcar-averaged lin e ratios with 
extinction corrections from the Schult heis et al.l (2009) map 



(dashed green line) and without any correction for extinction 
(dashed grey line). The [S III] 18.71 pm / fS ml 33.48 urn 
line ra tio is strongly affected by extinction, as ISimpson et al.l 
(2007) also found, and we therefore move the discussion of 
this line ratio to the Appendix. All of the other line ratios in 
Figure [10] are insensitive to the extinction correction, and will 
be discussed in the following sections. 

We illustrate the distribution of [Fe III] 22.93 pm I [Fe II] 
25.99 pm as a function of Galactic longitude in Figure [TOl 
This line ratio shows a broad peak centered at / ~ +0.15°, 
near the Quintuplet and Arches clusters of massive stars and 
the Radio Bubble. The [Fe III] 22.93 pm I [Fe II] 25.99 pm 
ratio is not sensitive to the hardness of the ionizing radia- 
tion field, as [Ne III] 15.56 pm I [Ne II] 12.81 ^m is (see 
below). Instead, it is sensitive to the ionization parameter, 
U, which is the r atio of the pho toionization rate to the re- 
combination rate (Contini 2009). This is because [Fe III] 
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22.93 /im, with an ionization potential of 16.2 eV, only 
arises in ionized gas, while [Fe II] 25.99 ^m, with an ion- 
ization poten tial of 7.9 eV, can ar ise both from neutral and 
ionized gas (Kauf man et alj 12006). Comparisons of mod- 
els with prev ious GC ISM observations have found —3 ^ 
Ior(U) < -1 dRodrfguez-Fernandez & Martfn-Pintadol 120051; 
lContinill2009h . 

4. COMPARISON WITH NEARBY GALAXIES 

In this section, we compare ratios of ionic lines mapped 
in the GC to ratios measure d in the Spitzer I nfrared Nearby 
Galaxies Survey (SINGS') (iDale et al.1 120091 The SINGS 
sample dKennicutt et alj 12003b contains 75 galaxies, evenly 
distributed among elliptical, spiral, and irregular galaxies, 
and including a range of nuclear activity (quiescent, starburst, 
LINER, Seyfert). The SINGS galaxies have a median distance 
of 9.5 Mpc, and the closest is at 0.6 Mpc. The IRS SH and LH 
slits correspond to 14 pc x33 pc and 32 pc x65 pc, respec- 
tively, for a galaxy at 0.6 Mpc; they cover 220 pc x520 pc 
and 510 pc x 1000 pc, respectively, for a galaxy at 9.5 Mpc. 
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By comparison, we bin our GC line data into 1 .5' x 1 .5' pixels 
(3.5 pc x3.5 pc). Our coadded spectrum of the CMZ is con- 
structed from spectra across a 210 pc x60 pc region, which is 
a good match to the SINGS spatial resolution. 

4. 1 . Radiation Field Hardness and Oxygen Abundance 

The [Ne III] 15.56 (im I [Ne II] 12.81 ^m line ratio is 
a useful indicator of the radiation field hardness for star- 
forming regions, as this ratio is higher when stars are hot- 
ter. The map ping results for this line ratio are included 
in Figure [101 As shown in this map, the highest exci- 
tation gas traced by [Ne III]/[Ne II] ratio is peaked in the 
Radio Bubble and Quintuplet cluster, where ionizing pho- 
tons from newly born massive stars in the Quintuplet clus- 
ter are most likely r esponsible for the hard radiation field. 
Simps on et al.l (120071) also found the same result along the / ~ 
0.1° stripe (see also Rodrigue z-Fernandez & Martm-Pintadol 
120051) . The mean line ratio (w 0.1-1) is consistent with 
recent bursts of massive star form ation in this region 
in the last few million years (e.g., iThornlev et all l2000t 
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Rodri guez-Fernandez & Martfn-Pin tado 2005), and is insen- 
sitive to a choice of foreground extinction corrections (see 
dashed lines in Figure ITOb. 

The [Ne III] / [Ne II] ratio also depends on the nebular oxy- 
ge n abundance, [ O/H], FigureQT]is a modification of Figure 4 
in iDale et alJ (|2009), which shows [Ne III]/[Ne II] measured 
from a sample of normal star-forming regions (both nuclear 
and extra-nuclear; open blue circles) in SINGS. As shown 
in Figure QT| the line ratio decreases as the nebular oxy- 
gen abundance increases in these normal extragalactic star- 
forming regions, because the UV stellar spectrum of ionizing 
hot stars depends on the (photo spheric) metal abundance. Ac- 
tive galactic nuclei (AGNs) are also displayed as red filled 
triangles, but they do not follow the line ratio vs. abundance 
trend observed among extra-galactic star-forming regions. 

The median [Ne III] / [Ne II] ratio (« 0.08) measured in 
the GC is shown as a filled diamond point in Figure QT| and 



the number distribution of GC spectra is displayed in the 
right panel. The vertical error bars indicate the interquartile 
range of line ratios measured in the GC (= 0.05-0.14). The 
line ratios from the coadded spectrum (Table |3), shown as 
an open box symbol, are consistent with the median value 
within the interquartile ranges. The [Ne III] / [Ne II] line 
ratio of the CMZ in Figure [TT] is plotted at the stellar oxy- 
gen a bundance, [12 + log (O/H)] = 9.04±0.19 (ICunha et all 
2007), which was derived from high-resolution infrared spec- 
troscopy of five luminous cool stars within 30 pc of the GC. 
Their abundance measurement is consistent with lDavies et al.1 
(2009), who measured [12 + l og(Q/H)1 = 9.09± 0.11 from a 
star that was not included in Cun ha et alJ (|2007). Figure [TT] 
demonstrates that the CMZ follows the abundance vs. [Ne III] 
/ [Ne II] trend observed in normal star-forming regions in 
nearby galaxies. This implies that the hardness of the stel- 
lar energy distribution in the GC is similar to those found in 
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FIG. 10.— Cont'd. 



nuclear and extra-nuclear star-forming environments, assum- 
ing that the nebular oxygen abundance [O/H] is close to what 
was measured from stars in the GC. 

4.2. AGN vs. Normal Star-Forming Activities in the GC 

Empirical evidence suggests that some mid-IR line ratios 
are useful diagnostic tools for di scriminating between normal 
star-forming regions an d AGNs dLutz et al.l fl998; Stur m et ail 
l2006tlDale et alJ2009l) . In this section, we use several of these 
emission-line diagnostics, such as [Si II] 34.82 /im / [S III] 
33.48 /im and [O IV] 25.89 ^m / [Ne II] 12.81 /im, to study 
properties of gas clouds in the CMZ. We inspect these line 
ratios in our GC spectra, both separately and as a whole, and 
compare results with those observed in other nearby galaxies. 

4.2.1. [Si II]/[S III] and [Fe II]/[Ne II] 

Figu re [T2l which is a modification of Figure 6 in lDale et al.l 
(2009), is one of such diagnostic tools, showing [Si II] 
34.82 fjm I [S III] 33.48 fim vs. [Fe II] 25.99 fim I [Ne II] 



12.81 /im from SINGS. Both nuclear and extra-nuclear star- 
forming regions are shown as blue circles, and AGNs are 
shown as red filled triangles. Previous theoretical work 
showed that the [Si II] 34.82 fim line can be emitted either 
from PDRs or H I I regions, while [S III ] 33.48 fim is mostly 
from H II regions (Kaufman et al. 2006). 

In Figure Q~2] normal star-forming regions show a strong 
correlation between [Si II] 34.82 fim I [S III] 33.48 fim and 
[Fe II] 25.99 fim I [Ne II] 12.81 fim, the latter being another 
useful line diagnostic of the ionization parameter. However, 
AGNs generally show higher values of these line ratios than 
values from normal star-forming nuclear and extra-nuclear re- 
gions in nearby galaxies, and a combination of these two IR 
line ratios separates AGNs and normal star-forming regions 
relatively well. Physical reasons for this empirical division 
are debated, but must include the low ionization potentials 
of Si + and Fe + ions (8.15 eV and 7.90 eV, respectively) and 
the higher ionization po tentials of Ne + an d S" 1 " 1 " (21 .56 eV and 
23.34 eV, respectively). Dale et ail (120091) proposed a number 
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of possible physical mechanisms for the observed separation 
between AGNs and normal star-forming regions, which in- 
clude (1) enhanced dust destruction and sublimation of refrac- 
tory elements (Si and Fe) in the harsh AGN environment, (2) 
extended low-ionization volumes in AGNs produced by X-ray 
photo-ionization processes, and/or (3) enhanced line emission 
from [Si II] 34.82 /im and [Fe II] 25.99 fim due to high gas 
density in PDRs and/or X-ray dominated regions in AGNs. 

Gray cross points in Figure Q~2]display line ratios measured 
from individual GC spectra, shown only for those detected at 
more than a 3er significance. Their number distributions are 
shown in a histogram on each axis. The filled diamond point 
represents the median line ratio from the individual spectra 
of the entire CMZ ([Si II] / [S III] = 1.21 and [Fe II] / [Ne II] 
= 0.06), with error bars indicating interquartile ranges. The 
values from the coadded CMZ spectrum are shown as an open 
box symbol; they are different from the median GC values be- 
cause we only plot 3er detections. Clearly, the vast majority 
of the GC measurements and the median GC values are found 
in an area that is mostly populated by normal star-forming 
regions in nearby galaxies. Line ratios on both axes are insen- 
sitive to extinction corrections. 

4.2.2. [O IV]/[Ne II] , [Fe II]/[0 IV], and [Ne IlI]/[Ne II] 

Another set of mid-IR line diagnostics to separate AGNs 
from normal star-forming regions are displayed in Figure [13] 
which plots [O IV] 25.89 /im / [Ne II] 12.81 [im vs. [Fe II] 
25.99 /im / [O IV] 25.89 /mi, and in Figure[l4] which shows 
[O IV] 25.89 /im I [Ne II] 12.81 fim vs. [Ne III] 15.56 fim I 
[Ne II] 12.81 fim. Both diagrams, which involve line emis- 
sion from [O IV] 25.89 /im, separate AGNs (red filled tri- 
angles) from normal star-forming re gions (blue open circles) 
relatively well; see iDale et all (120091) for the SINGS targets 



used in these diagrams. In addition, two diagonal dashed 
lines in Figure [13] show approximate divisions between star- 
bursts, LINERs (Low Ionization Nuclear Emission Regions') 
and pure shocks, and Seyfert galaxies from S turm et al.l (2006. 
see their Figure 2). Figure [14] shows that there are a few 
normal star-forming regions from SINGS (blue circles) with 
[O IV] 25.89 [an I [Ne II] 12.81 fim greater than - 0.10. 
These are low-metallicity star-forming regions, as revealed 
from the comparison between Figure Q~T]and Figure [14] where 
strong [O IV] 25.89 fim emission is due to a harder spec- 
trum for the stellar ionizing photons in low-metallicity envi- 
ronments. 

Gray cross points in Figures [T3TjT4l represent our measure- 
ments from individual GC spectra. As in Figure Q~2] we only 
included emission lines with more than 3a detections. Num- 
ber distributions of ionic line ratios are shown on each axis. 
The median values from the individual spectra of the CMZ 
are indicated by a grey diamond point; error bars represent 
the interquartile ranges of distributions of individual GC spec- 
tra. The line ratios for the coadded GC spectrum is displayed 
as an open box. Figures [1311141 show that the emission-line 
properties of the CMZ are similar to those of extragalactic 
star-forming regions (blue circles) and starburst galaxies (up- 
per left corner in Figure [13] fenced with a dashed diagonal 
line). This agrees with our conclusion based on [Si II] / [S III] 
vs. [Fe II] / [Ne II] diagnostics in Figure [T2l 

The line ratios for the coadded spectrum of the CMZ in 
Figure [H([0 IV]/[Ne II]- 0.01 and [Ne III]/[Ne n]~ 0.08) 
are found within the range observed a mong starburst ga laxies 
jLutzetalJU998L see their Figure 2). iLutz et al.1 (1998) sug- 
gested that shocks and/or hot stars are the likely origins of 
JO IV] observed in th ese galaxies, rather than buried AGNs. 
Simps on" et al.l (2007) also reported widespread detection of 
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[O IV] i n the GC; they attributed this to shoc ks in the CMZ 
(see also IContinlll2009l) . ISimpson"et al. (2007) further argued 
for the existence of shocked gas inside of the Radio Bubble 
based on the high Fe abundance there, since grain destruction 
in shocks can return Fe to the ISM. 

If the above result implies the importance of shocked gas in 
the GC, widespread detections of highly ionized ions such as 
[O IV] (Figure [7]i may support a large scale origin of turbu- 
lent motions in the GC. In other words, the relatively uniform 
intensity distributions of these lines in our mapping suggest 
that conditions for producing the lines are common through- 
out the CMZ. One likely explanation is a dissipation of su- 
personic turbulence (> 10 km/s in the CMZ) induced by dif- 
ferential Gal actic rotation or shearing motion of cl ouds in 
the GC (e.g.. iWilson eTaTlll982t iGuesten et al.l[l985l) . How- 
ever, shocks created by massive stars are also likely; see 
Rodriguez-Fernande z et al.1 (|2004) for discussion of other po- 
tential heating mechanisms, such as cosmic rays, magnetic 
heatings, or X-ray dominated regions. 

While most of the GC points (gray crosses) in Figure [T21 are 
found in an area occupied by normal star-forming regions, in- 
cluding the coadded spectrum of the CMZ, about 10% of the 
detections fall into AGN territory. As shown in the mapping 
results in Figure [10] both [Si II] 34.82 /im I [S III] 33.48 /zm 
and [Fe II] 25.99 /im / [Ne II] 12.81 /im line ratios are peaked 
near Sgr A and the north-western rim of Sgr B. The scatter 
of the [Fe II] 25.99 fim I [Ne II] 12.81 /im points in the mov- 
ing average plot is larger than that of [Si II] 34.82 /im / [S III] 
33.48 /im because [Fe II] 25.99 /im is relatively weak and 
blended with [O IV] 25.89 /im. Nevertheless, the observed 
systematic trend in Galactic longitude (right panels in Fig- 
ure [Tot is larger than the random scatter of data points, which 
means that those GC points in the upper right corner in Fig- 
urell2lare not entirely produced by a statistical fluctuation. 

Figures [T3l and fl4l show that about 10% of individual GC 
spectra fall in the AGN (red filled triangles) region of these 
line ratio diagrams, almost independent of the choice of ex- 
tinction correction. Note again that a few SINGS targets from 
extragalactic star-forming regions are also found in the area 
denoted as AGNs in Figure [14] because of their low metal 
abundance (see above). On the other hand, the GC spectra 
with strong [O IV] 25.89 /im / [Ne II] 12.81 /zm are not due to 
a low metallicity environment, because they follow the AGN 
sequence in Figure [14] rather than that of metal-poor, star- 
forming reg ions, and because the GC has a supers olar oxygen 
abundance (ICunha et al.ll2007l: iDavies et alJl2009h . 

If high values of [O IV] 25.89 /zm / [Ne II] 12.81 /im and 
low values of [Fe II] 25.99 /im / [O IV] 25.89 /im in Figures[T3l 
are due to ionization by a power-law continuum source (i.e., 
an AGN), then the empirical division in Figure Q~2] leads to 
a conclusion that these same spectra should also have high 
values of the ratios of [Si II] 34.82 /jm / [S III] 33.48 /zm and 
[Fe II] 25.99 /im / [Ne II] 12.81 /zm. To answer the question 
of whether AGN-like points in Figures [T2T - [l4l are from the 
same individual GC spectra, we plot [Si II] 34.82 /zm / [S III] 
33.48 /zm vs. [OIV] 25.89 /im/ [Nell] 12.81 /tm in Figure[B] 
We also plot these line ratios measured in the coadded CMZ 
spectrum, and compare to th e line ratios m easured in star- 
forming galaxies and AGNs (Dal e et aLll2009l) . 

We find that only 3% of GC points fall in the AGN-like 
region of Figure [15] Virtually all GC data points have ratios 
of [Si II] 34.82 /im / [S III] 33.48 /zm and [O IV] 25.89 /zm 
/ [Ne II] 12.81 /zm agreeing with the line ratios observed in 
star-forming galaxies. The few points that are not similar to 



star-forming galaxies are those at high Galactic latitudes (/ ~ 
+0.4°, b ~ +0.2° and / - +0.6°, b ~ +0.2°). 

We considered whether the positions with AGN-like line 
ratios could be due to having been irradiated by Sgr A* 
if it was in a more active state in the past. This idea 
has been suggested to explain bright Fe I Ka emission at 
6.4 keV in Sgr B2 and other molecular c l ouds in the GC 
dSunyaev et all 11993b iKovama et ail 119961: iMurakami eTai] 
1200 ID . The Fe I Ka emission in the GC varies spatially 
and temporily (iTerrier et al.l 120101) . ICapelli et all (120121) an- 
alyzed the Fe I Ka emission to derive a light curve for 
the X-ray luminosity of Sgr A* over the last several hun- 
dred years. Apparent superluminal motio n has also bee n 
observed in 6.4 keV Fe I Ka emission dPonti etalJUoiOh . 
Others argue, however, that widespread F e I Ka emission in 
the GC can be explained by co smic r avs (lYu sef-Zad eh et al.1 
[20071 120T1 IChernyshov et al.ll2012t Ifatischeff et al.1 120121) . 
lYusef-Zadeh et al] (120 1 3h show, in their Figure 8a, an im- 
age of the equivalent width of Fe I Ka emission, covering 
-0.8° < I < +0.7° and ±0.4° in b. We compared this image 
to ourFigure[lO] an image of the [O IV]/[Ne II] line ratio in the 
GC ISM. We find no correlation between GC ISM positions 
with AGN-like line ratios and regions of enhanced Fe I Ka 
emission. Our data do not support - but neither do they rule 
out - the idea that strong 6.4 keV Fe I Ka emission is caused 
by higher X-ray luminosity of Sgr A* in the past. 

We conclude that our observations of mid-IR line emis- 
sion in the central 210 pc x60 pc of the Galaxy show no 
evidence of excitation by a power-law continuum. The few 
points that are not similar to star-forming galaxies are out- 
liers, such as those excited by very hot sources (e.g., Wolf- 
Rayet stars, planetary nebulae, or X-ray binaries), rather than 
be longing to an AGN-lik e tre nd. Our results agree with those 
of lSimpson et all (120071) and IContinil (120091) . who conclude 
that the mid-IR line ratios in a smaller Spitzer data set can 
be explained by photoionization by hot stars, combined with 
shocks. 

Figure [15] shows a significant offset between the median 
value of [O IV] 25.89 /zm / [Ne II] 12.81 /zm measured from 
individual GC spectra and the value from the coadded CMZ 
spectrum. Emission from [Ne II] 12.81 /zm, [Si II] 34.82 /zm, 
and [S III] 33.48 /tm is detected in virtually all individual GC 
spectra. The offset, then, is due to faint [O IV] 25.89 /im 
emission which is missing from the number distribution for 
individual spectra (we only plot 3er detections) but which 
contributes to the coadded spectrum. Because the [O IV] 
25.89 /im/ [Ne II] 12.81 /zm line ratio from the coadded spec- 
trum is at the lower end of the range of GC line ratios on Fig- 
ure Q3] we can be confident that excluding regions of faint 
[O IV] 25.89 /zm emission does not change our conclusion 
that the GC is similar to star-forming galaxies and not similar 
to AGNs. 

5. SUMMARY 

We present a mid-IR spectroscopic survey of 428 posi- 
tions in the ISM of the CMZ using the SpitzerflRS, and con- 
struct mid-IR emission line maps for several forbidden and 
molecular hydrogen lines over the CMZ. We derive both line 
strengths and radial velocities from individual lines, and com- 
pute line flux ratios as a probe of physical conditions in the 
GC. Our mapping is superior to previous survey results in 
terms of the total area covered in the GC. We also construct a 
CMZ spectrum by coadding individual spectra after correct- 
ing each for extinction. 
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Mid-IR emission lines from the pure rotational transitions 
of molecular hydrogen, S(0), S(l), and S(2), are observed in 
almost all lines of sight to the CMZ. Their intensity distribu- 
tion in the GC is relatively uniform; their radial velocity dis- 
tributions are poorly correlated with those from ionic species, 
with worse correlation for S(0) than S(2). We view this as evi- 
dence that most of the H2 S(0) emission, and some of the S(l) 
and S(2) emission, arises from PDRs along the line of sight to 
the GC, rather than PDRs associated with ionized gas in the 
GC. 

The radiation field hardness traced by the [Ne III] 15.56 [im 
I [Ne II] 12.81 fim line ratio indicates that the highest exci- 
tation gas clouds are found in the Radio Bubble region and 
Quintuplet cluster, and the mean value in the GC is consistent 
with a recent burst of star formation in the last few million 
years. The hardness of the ionization spectrum from hot stars 
is tied to the metal abundance of the star-forming clouds, and 
our GC spectra, combined with the published GC stellar oxy- 
gen abundance, show that the hardness of the GC exciting ra- 
diation is similar to that found in normal star-forming regions 
in nearby galaxies. 

We present mid-IR line-ratio diagrams such as [Si II] 
34.82 /im / [S III] 33.48 fj,m vs. [Fe II] 25.99 fim I [Ne II] 
12.81 /im, [O IV] 25.89 fim / [Ne II] 12.81 /jm vs. [Fe II] 
25.99 fim I [O IV] 25.89 /xm, [O IV] 25.89 /im / [Ne II] 
12.81 /Lim vs. [Ne III] 15.56 fxm I [Ne II] 12.81 /im, and [Si II] 



34.82 /im / [S III] 33.48 /im vs. [O IV] 25.89 /mi I [Ne II] 
12.81 /im. We compare properties of individual GC spectra 
to those observed in nearby extragalactic star-forming regions 
and AGNs. These diagrams show that the individual GC spec- 
tra and the mean GC spectrum are consistent with normal star 
forming activity, where emission from highly ionized species 
such as [O IV] is likely produced by shocks and/or turbulence 
prevalent in the CMZ clouds. Our GC line ratios do not agree 
with line ratios observed for LINER galaxies or AGNs. 
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APPENDIX 
[S III] 18.71 fiM I [S III] 33.48 /iM 

Line ratios from the same ionized species such as [S III] 18.71 11 ml 33.48 /mi are generally sensitive to the electron gas density 
(n e ) with a weak dependence on the electron gas temperature (e.g., Rubin 1989; Martin-Hernandez et al. 2002; Dale et al. 2009). 
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But [S III] 18.71 fj,m is found in the middle of the strong and broad 18 /im silicate feature, while JS III] 33.48 fim is almost 
unaffected by foreground absorption (see the bottom panel in Figure|2]i. Thus we, like Sim pson et al.l (120071) . find that the [S III] 
18.71 /im / 33.48 /im ratio is most useful for checking our foreground extinction correction. 

The grey dashed line for the [S III] 18.71 /im and 33.48 /im ratio in the right panel of Figure [TOlindicates the observed line 
intensity ratios without any extinction correction. These values are located below the lower limit set by the theoretical value 
(0.4; dotted horizontal li ne) as n e — > cm" " 3 (Rubin 1989)B On the other hand , mean line ratios after extinction corrections using 
either the Simp son et al.1 (120071) method (red solid line) or the lSchultheis et ai] (2009) map (green dashed line) are higher than the 
theoretical limit, clearly showing a strong dependence of the line ratio on the foreground extinction correction. 

Over most parts of the GC, the [S III] 18.71 /im / [S III] 33.48 fim ratios are close to unity after the extinction correction. This 
line ratio value corresponds to n e ~ 500 cm -3 dMartm-Hemand ez et al.ll2002t iDale et al.ll2009l) . Mean line ratios from either of 
the extinction estimation techniques used above are in excellent agreement with each other, and are relatively uniform across the 
CMZ. However, these average line ratios become lower than the minimum theoretical value (w 0.4) near Sgr B2. Sgr B2 has 
very high extinction (Figure B), so we view this as due to unrecognized systematic errors in our extinction e stimates from the 
ISimpson et al.l (120071) and Schulth eis et al.l J2009) approaches that greatly underestimated Tgj in this region. iGoicoechea et al.1 
d2004l) found n e ~ 240 cm 1 in Sgr B2 based on ISO observations. 



5 New cross sections of S III (Hudson et alj EOTJ) produce a lower limit of ^ Q 5 for ^ ^ ratio (J simpsQn 1Ql ^ private communication) . 



